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^ ■ 1 Historical Summary 

Cosmological dark matter in the form of neutrinos with masses of up to a 
few electron volts is known as hot dark matter. In 1979-83, this appeared 
■ to be perhaps the most plausible dark matter candidate. Such HDM models 

of cosmological structure formation led to a top-down formation scenario, in 
which superclusters of galaxies are the first objects to form, with galaxies 
and clusters forming through a process of fragmentation. Such models were 
abandoned when it was realized that if galaxies form sufficiently early to agree 
I with observations, their distribution would be much more inhomogeneous 

than it is observed to be. Since 1984, the most successful structure formation 
I models have been those in which most of the mass in the universe is in the 

i form of cold dark matter (CDM). But mixed models with both cold and hot 

I dark matter (CHDM) were also proposed in 1984, although not investigated 

in detail until the early 1990s. 

The recent atmospheric neutrino data from Super-Kamiokande provide 
strong evidence of neutrino oscillations and therefore of non-zero neutrino 
O I mass. These data imply a lower limit on the HDM (i.e., light neutrino) con- 

tribution to the cosmological density il^, ^ 0.001 — almost as much as that 
^ I of all the stars in the centers of galaxies — and permit higher il^. The "stan- 

dard" COBE-normalized critical-matter-density (i.e., f2m = 1) CDM model 
has too much power on small scales. It was discovered in 1992-95 that CDM 
i with the addition of neutrinos with total mass of about 5 eV, corresponding 

I to ~ 0.2, results in a much improved fit to data on the nearby galaxy and 

cluster distribution. Indeed, the resulting Cold -I- Hot Dark Matter (CHDM) 
cosmological model is arguably the most successful i7„, = 1 model for struc- 
ture formation ||l|,|,|,|. 

However, other recent data have begun to make a convincing case for 
0.3 ^ ilrn ~ 0.5. In light of all these new data, several authors have considered 
whether cosmology still provides evidence favoring neutrino mass of a few eV 
in flat models with cosmological constant Ha = 1 — f^m- The conclusion 
is that the possible improvement of the \ow-ilm flat (ylCDM) cosmological 
models with the addition of light neutrinos appears to be rather limited, 
but that ylCHDM models with ^0.1 may be consistent with presently 
available data. Data expected soon may permit detection of such a hot dark 
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matter contribution, or alternatively provide stronger upper limits on J?^ and 
neutrino masses. 

2 Hot, Warm, and Cold Dark Matter 

Hot DM refers to particles, such as neutrinos, that were moving at nearly the 
speed of light at redshift z ~ 10^ (or time t ~ 1 yr), when the temperature 
T ~ 3 X 10^ eV and the cosmic horizon first encompassed lO^^Mo, the amount 
of dark matter contained in the halo of a large galaxy like the Milky Way. 
Hot DM particles must also be still in thermal equilibrium after the last 
phase transition in the hot early universe, the QCD confinement transition, 
which presumably took place at Tqcd ~ 10^ MeV. Hot DM particles have 
a cosmological number density roughly comparable to that of the microwave 
background photons, which (as we will see shortly) implies an upper bound 
to their mass of a few tens of eV. This then implies that free streaming of 
these relativistic particles destroys any fluctuations smaller than supercluster 
size, ~ IO^^Mq. 

The "hot," "warm," "cold" DM terminology was introduced in 1983 §,§. 
Warm DM particles interact much more weakly than neutrinos. They decou- 
ple (i.e., their mean free path first exceeds the horizon size) at T ^ TqcDi 
and are not heated by the subsequent annihilation of hadronic species. Con- 
sequently their number density is roughly an order of magnitude lower, and 
their mass an order of magnitude higher, than hot DM particles. Fluctua- 
tions corresponding to sufficiently large galaxy halos, ^ lO^^M©, could then 
survive free streaming. In theories of local supersymmetry broken at ^ 10^ 
GeV, gravitinos could be DM of the warm variety (7|J|,^. Other warm dark 
matter candidates are also possible, of course, such as right-handed neutri- 
nos Warm DM does not fit the observations if f2,„ = 1 but for 
low fl^ some have suggested that it may be worth reconsidering, to avoid 
some possible problems of Cold DM ||l^,|l3|. However, the cutoff in the power 
spectrum P(fc) at large fc implied by WDM will also inhibit the formation 
of small dark matter halos at high redshift. But such small halos are pre- 
sumably where the first stars form, which produce metals rather uniformly 
throughout the early universe as indicated by observations of the Lyman a 
forest (neutral hydrogen clouds seen in absorption in quasar spectra). 

Cold DM consists of particles for which free streaming is of no cosmo- 
logical importance. Two different sorts of cold DM consisting of elementary 
particles have been proposed, heavy thermal remnants of annihilation such 
as supersymmetric neutralinos, and a cold Bose condensate such as axions. 
A universe where the matter is mostly cold DM and there is a large cosmo- 
logical constant looks very much like the one astronomers actually observe, 
and this low-J7™ ylCDM model llitl is the current favorite model for structure 




Hot Dark Matter 3 



3 Galaxy Formation with Hot DM 

The standard hot DM candidate is massive neutrinos , ahhough 

other, more exotic theoretical possibilities have been suggested, such as a 
"majoron" of nonzero mass which is lighter than the lightest neutrino species, 
and into which all neutrinos decay. Neutrinos appeared to be an attractive 
DM candidate because of the measurement of an electron neutrino mass of 
about 30 eV in 1980 fH]. This coincided with the improving CMB Hmits 
on the primordial fluctuation amplitude, which forced Zel'dovich and other 
theorists to abandon the idea that all the dark matter could be made of 
ordinary baryonic matter. The version of HDM that they worked out in detail, 
with adiabatic Gaussian primordial fluctuations, became the prototype for 
the subsequent i?,,, = 1 CDM theory. 

3.1 Mass Constraints 

Direct measurements of neutrino masses have given only upper limits (see also 
the chapter by Robertson and Wilkerson). A secure upper limit on the elec- 
tron neutrino mass is roughly 15 eV. The Particle Data Group ||2^ notes that 
a more precise limit cannot be given since unexplained effects have resulted 
in significantly negative measurements of m{i'e)^ in tritium beta decay ex- 
periments. However, this problem is at least partially resolved, and the latest 
experimental upper limits on the electron neutrino mass are 2.8 eV from the 
Mainz and 2.5 eV from the Troitsk tritium beta decay experiments 
(both 95% C.L.). There is an upper limit on an effective Majorana neutrino 
mass of ~ 1 eV from neutrinoless double beta decay experiments The 
upper limits from accelerator experiments on the masses of the other neutri- 
nos are m{iy^,) < 0.17 MeV (90% CL) and m(i/^) < 18 MeV (95% CL) |2|,|3, 
but since stable neutrinos with such large masses would certainly "overdose 
the universe" (i.e., contribute such a large cosmological density that the uni- 
verse could never have attained its present age), cosmology implies a much 
lower upper limit on these neutrino masses. 

Before going further, it will be necessary to discuss the thermal history 
of neutrinos in the standard hot big bang cosmology in order to derive the 
corresponding constraints on their mass. Left-handed neutrinos of mass < 1 
MeV remain in thermal equilibrium until the temperature drops to Ti^d, at 
which point their mean free path first exceeds the horizon size and they es- 
sentially cease interacting thereafter, except gravitationally p8| . Their mean 
free path is, in natural units {h ^ c ^ 1), Ai, [cr^riei]"^ ~ [{GpT'^){T^)]^''- , 
where Gp « 10^*^ GeV~^ is the Fermi constant that measures the strength of 
the weak interactions. The horizon size is \h {Gp)~^^'^ ^ MpiT~^, where 
the Planck mass Mpg = G'^/^ = 1.22 x lO^^ GeV. Thus Xh/X^ - {T/T^df, 
with the neutrino decoupling temperature 

T,d Mpl'^'Gp^''' ~ 1 MeV . (1) 
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After T drops below i MeV, e+e~ annihilation ceases to be balanced by 
pair creation, and the entropy of the e^e~ pairs heats the photons. Above 
1 MeV, the number density n^i of each left-handed neutrino species and its 
right-handed antiparticle is equal to that of the photons, n^, times the factor 
3/4 from Fermi versus Bose statistics. But then e^e~ annihilation increases 
the photon number density relative to that of the neutrinos by a factor of 
11/4. ([]) As a result, the neutrino temperature T^,o = {'^/llY^'^T^.o- Thus 
today, for each species, 

riu^o = \-^ "7,0 = 109 6'3cm-3 , (2) 

where 6 = (To/2.7K). With the cosmic background radiation temperature 
To = 2.728 ± 0.004 K measured by the FIRAS instrument on the COBE 
satellite [||, r^,o = 1-947 K and n^,o = 112 cm"^ 
Since the present cosmological matter density is 

p,„ = Qpc = 10.54 f2„/i2keVcm~3 , (3) 

it follows that 

m„ < prrjn^.o < 96 Qmh^e^^eY « 93 Hmh^eY , (4) 

i 

where the sum runs over all neutrino species with M^i < 1 MeV. (Heavier 
neutrinos will be discussed in the next paragraph.) Observational data imply 
that Hmh^ « 0.1 - 0.3, since fim ~ 0.3 - 0.5 and h « 0.65 ± 0.1 |l|l. Thus 
if all the dark matter were light neutrinos, the sum of their masses would be 
« 9 - 28 eV. 

In deriving eq. (jj), we have been assuming that all the neutrino species 
are light enough to still be relativistic at decoupling, i.e. lighter than an MeV. 
The bound shows that they must then be much lighter than that. In the 
alternative case that a neutrino species is nonrclativistic at decoupling, it 
has been shown |^ , ^l|j3^ , |33| , ^ that its mass must then exceed several GeV, 
which is not true of the known neutrinos (i/g, i/^, and Vr)- (One might at 
first think that the Boltzmann factor would sufhciently suppress the number 

^ In the argument giving the 11/4 factor, the key ingredient is that the entropy 
in interacting particles in a comoving volume Si is conserved during ordinary 
Hubble expansion, even during a process such as electron-positron annihilation, 
so long as it occurs in equilibrium. That is. Si = gi{T)N-f{T) — constant, where 
A''^ = n-yV is the number of photons in a given comoving volume V, and gi = 
[gs + f fff)/ is the effective number of helicity states in interacting particles (with 
the factor of | reflecting the difference in energy density for fermions versus 
bosons). Just above the temperature of electron-positron annihilation, gi = g-^i + 
|x(7e = 2-f|x4=^; while below it, gi ~ g-i = 2. Thus, as a result of 
the entropy of the electrons and positrons being dumped into the photon gas at 
annihilation, the photon number density is thereafter increased relative to that 
of the neutrinos by a factor of 11/4. 
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density of neutrinos weighing a few tens of MeV to allow compatibility with 
the present density of the universe. It is the fact that they "freeze out" of 
equilibrium well before the temperature drops to their mass that leads to 
the higher mass limit.) We have also been assuming that the neutrino chem- 
ical potential is negligible, i.e. that — Tipj <C n^- This is very plausible, 
since the net baryon number density (rib — ni) H 10~^n^, and big bang nu- 
cleosynthesis restricts the allowed parameters |^ (see also the chapter by 
Fuller). 



3.2 Phase Space Constraint 

We have just seen that light neutrinos must satisfy an upper bound on the 
sum of their masses. But now we will discuss a lower bound on neutrino mass 
that arises because they must be rather massive to form the dark matter 
in galaxies, since their phase space density is limited by the Pauli exclusion 
principle. A slightly stronger bound follows from the fact that they were not 
degenerate in the early universe. 

The phase space constraint ||3^] follows from Jeans 's theorem in classical 
mechanics to the effect that the maximum 6-dimensional phase space density 
cannot increase as a system of collisionless particles evolves. At early times, 
before density inhomogenitites become nonlinear, the neutrino phase space 
density is given by the Fermi-Dirac distribution 



pc 



(5) 



where here h is Planck's constant and = 2 for each species of left-handed 
u plus right-handed D. Since momentum and temperature both scale as red- 
shift z as the universe expands, this distribution remains valid after neu- 
trinos drop out of thermal equilibrium at ^ 1 MeV, and even into the 
nonrelativistic regime < [ p8| . The standard version of the phase 
space constraint follows from demanding that the central phase space density 
9[2(27r)^/^Gr^CTm^]~^ of the DM halo, assumed to be an isothermal sphere 
of core radius Tc and one-dimensional velocity dispersion a, not exceed the 
maximum value of the initial phase space density n^(0) — g^/2h^. The result 
is 

m.>(120eV)(^^^j (^j (y) ■ 

The strongest lower limits on riii, follow from applying this to the small- 
est galaxies. Both theoretical arguments regarding the dwarf spheroidal (dS) 
satellite galaxies of the Milky Way and data on Draco, Carina, and Ursa 
Minor made it clear some time ago that dark matter dominates the gravi- 
tational potential of these dS galaxies, and the case has only strengthened 
with time p8| . The phase space constraint then sets a lower limit > 
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500 eV, which is completely incompatible with the cosmological constraint 
eq. (^). However, this argument only excludes neutrinos as the DM in cer- 
tain small galaxies; it remains possible that the DM in these galaxies is (say) 
baryonic, while that in larger galaxies such as our own is (at least partly) 
light neutrinos. A more conservative phase space constraint was obtained for 
the Draco and Ursa Minor dwarf spheroidals |Q , but the authors concluded 
that neutrinos consistent with the cosmological upper bound on nii, cannot 
be the DM in these galaxies. A similar analysis applied to the gas-rich low- 
rotation- velocity dwarf irregular galaxy DDO 154 [4l| gave a limit rrii, > 94 
eV, again inconsistent with the cosmological upper bound. 



3.3 Free Streaming 

The most salient feature of hot DM is the erasure of small fluctuations by 
free streaming. Thus even collisionless particles effectively exhibit a Jeans 
mass. It is easy to see that the minimum mass of a surviving fluctuation 
is of order Mp^/ml ji^ , ^ . Let us suppose that some process in the very 
early universe — for example, thermal fluctuations subsequently vastly in- 
flated in the inflationary scenario — gave rise to adiabatic fluctuations on all 
scales. In adiabatic fluctuations, all the components — radiation and matter 
— fluctuate together. Neutrinos of nonzero mass stream relativistically 
from decoupling until the temperature drops to T m^, during which time 
they traverse a distance d,y = Rh{T = m,y) ~ Mpim~^. In order to survive 
this free streaming, a neutrino fluctuation must be larger in linear dimension 
than dj/. Correspondingly, the minimum mass in neutrinos of a surviving 
fluctuation is Mj^i, ~ dlm,jn„{T = m„) ~ dlmf, ~ Mp^m'^. By analogy 
with Jeans's calculation of the minimum mass of an ordinary fluid pertur- 
bation for which gravity can overcome pressure, this is referred to as the 
(free-streaming) Jeans mass. 

A more careful calculation gives 

= 4:l{mj30eY)-\l + z)-^Mpc , (7) 

that is, d^ = 41(m^/30eV)~^ Mpc in comoving coordinates, and correspond- 
ingly 

Afj„ = 1.77 M|,^m;2 = 3.2 x 10^^(m^/30 eV)-^M0 , (8) 

which is the mass scale of superclusters. Objects of this size are the first to 
form in a z^-dominated universe, and smaller scale structures such as galaxies 
can form only after the initial collapse of supercluster-size fluctuations. 

When a fluctuation of total mass ~ lO^^M© enters the horizon at z ~ 10**, 
the density contrast Jrb of the radiation plus baryons ceases growing and 
instead starts oscillating as an acoustic wave, while that of the massive neu- 
trinos Si, continues to grow linearly with the scale factor R= {1 + z)~^ since 
the Compton drag that prevents growth of (Jrb does not affect the neutrinos. 
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By recombination, at ^ 10^, ^rb/'^,/ ^ 10"^, with possible additional sup- 
pression of (5rb by Silk damping. Thus the hot DM scheme with adiabatic 
primordial fluctuations predicts small-angle fluctuations in the microwave 
background radiation that are lower than in the adiabatic baryonic cosmol- 
ogy, which was one of the reasons HDM appealed to Zel'dovich and other 
theorists. Similar considerations apply in the warm and cold DM schemes. 
However, as we will discuss in a moment, the HDM top-down sequence of 
cosmogony is wrong, and with the COBE normalization hardly any structure 
would form by the present. 

In numerical simulations of dissipationless gravitational clustering start- 
ing with a fluctuation spectrum appropriately peaked at A ~ d,y (reflecting 
damping by free streaming below that size and less time for growth of the 
fluctuation amplitude above it), the regions of high density form a network of 
filaments, with the highest densities occurring at the intersections and with 
voids in between . The similarity of these features to those seen 

in observations was cited as evidence in favor of HDM Bfll . 



3.4 Problems writh v DM 

A number of potential problems with the neutrino dominated universe had 
emerged by about 1983, however. 

• From studies both of nonlinear clusterin g p^JSOt (comoving length scale 
A 10 Mpc) and of streaming velocities [^1[ in the linear regime (A > 10 
Mpc), it follows that supercluster collapse must have occurred recently: 
^sc — 0-5 is indicated and in any case z^^ < 2 However, the best 
limits on galaxy ages coming from globular clusters and other stellar 
populations indicated that galaxy formation took place before z « 3. 
Moreover, if quasars are associated with galaxies, as is suggested by the 
detection of galactic luminosity around nearby quasars and the apparent 
association of more distant quasars with galaxy clusters, the abundance 
of quasars at z > 2 was also inconsistent with the "top-down" neutrino 
dominated scheme in which superclusters form first: Zsc > Zgaiaxies- 

• Numerical simulations of the nonlinear "pancake" collapse taking into 
account dissipation of the baryonic matter showed that at least 85% of 
the baryons are so heated by the associated shock that they remain unable 
to condense, attract neutrino halos, and eventually form galaxies [|5|,|5^. 
This was a problem for the hot DM scheme for two reasons. With the 
primordial nucleosynthesis constraint Qb ^ 0.1, there would be difficulty 
having enough baryonic matter condense to form the luminosity that 
we actually observe. And, where are the X-rays from the shock-heated 
pancakes [|3)? 

• The neutrino picture predicts |Q that there should be a factor of ~ 5 
increase in M/Mi, between large galaxies (Af ~ IO^^Mq) and large clus- 
ters (A/ > lO^'^A^o), since the larger clusters, with their higher escape 
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velocities, are able to trap a considerably larger fraction of the neutrinos. 
Although there is some indication that the mass-to-light ratio M/L in- 
creases with Af, the ratio of total to luminous mass M/Mi^^ is probably 
a better indicator of the value of M/Mi,, and it is roughly the same for 
galaxies with large halos and for rich clusters. 

These problems, while serious, would perhaps not have been fatal for the 
hot DM scheme. But an even more serious problem for HDM arose from 
the low amplitude of the CMB fluctuations detected by the COBE satellite, 
(AT/T)rms = (1-1 ± 0.2) X 10~^ smoothed on an angular scale of about 10° 
[ p5| . Although HDM and CDM both have the Zel'dovich spectrum shape 
(P(fc) oc k) in the long-wavelength limit, because of the free-streaming cutoff 
the amplitude of the HDM spectrum must be considerably higher in order to 
form any structure by the present. With the COBE normalization, the HDM 
spectrum is only beginning to reach nonlinearity at the present epoch. 

Thus the evidence against standard hot DM is convincing. At very least, 
it indicates that structure formation in a neutrino-dominated universe must 
be rather more complicated than in the standard inflationary picture. 

The main alternative that has been considered is cosmic strings plus 
hot dark matter. Because the strings would continue to seed structure up 
until the present, and because these seeds are in the nature of rather lo- 
calized fluctuations, hot DM would probably work better with string seeds 
than cold DM. However, strings and other cosmic defect models are now es- 
sentially ruled out [ p6|j57[ | because they predict that the cosmic microwave 
background would have an angular power spectrum without the pronounced 
(doppler/acoustic/Sakharov) peak at angular wavenumber I ~ 220 that now 
appears to be clearly indicated by the data, along with secondary peaks at 
higher I. 

4 Cold plus Hot Dark Matter and Structure 
Formation: = 1 

Even if most of the dark matter is of the cold variety, a little hot dark matter 
can have a dramatic effect on the predicted distribution of galaxies. In the 
early universe, the free streaming of the fast-moving neutrinos washes out 
any inhomogeneities in their spatial distribution on the scales that will later 
become galaxies. If these neutrinos are a significant fraction of the total mass 
of the universe, then although the density inhomogeneities will be preserved 
in the cold dark matter, their growth rates will be slowed. As a result, the 
amplitude of the galaxy-scale inhomogeneities today is less with a little hot 
dark matter than if the dark matter is only cold. (With the tilt n of the 
primordial spectrum Pp{k) = Ak" fixed — which as we discuss below is not 
necessarily reasonable — the fractional reduction in the power on small scales 
is AP/P « 8f2y/ Qm |5^. See Fig. 1 for examples of how the power spectrum 
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P{k) is affected by the addition of hot dark matter in J?^ = 0.4 flat cos- 
mologies.) Since the main problem with f2„ = 1 cosmologies containing only 
cold dark matter is that the amplitude of the galaxy-scale inhomogeneities 
is too large compared to those on larger scales, the presence of a little hot 
dark matter appeared to be possibly just what was needed. And, as was men- 
tioned at the outset, a CHDM model with Sim = 1, f2,y — 0.2, and Hubble 
parameter h = 0.5 is perhaps the best fit to the galaxy distribution in the 
nearby universe of any cosmological model. The effects of the relatively small 
amount of hot dark matter in a CHDM model on the distribution of matter 
compared to a purely CDM model arc shown graphically in |^ ; cf . also |^ . 
As expected, within galaxy halos the distribution of cold and hot particles 
is similar. But the hot particles are more widely distributed on larger scales, 
and the hot/cold ratio is significantly enhanced in low-density regions. 

The first step in working out the theory of structure formation is to use 
linear perturbation theory, which is valid since cosmic microwave background 
measurements show that density fluctuations are small at the redshift of 
recombination, Zr ^ 10'^. The most extensive early calculations of this sort 
were carried out by Holtzman [l6l|j6^ , who concluded that the most promising 
cosmological models were CHDM and ylCDM . The most efficient method 
of computing the linear evolution of fluctuations now is that used in the 
CMBFAST code 116^ . An alternative Monte Carlo treatment of the evolution 



of neutrino density fluctuations was given by 1 65 , but the differences from the 
usual treatment appear to be small. Detailed analytic results have been given 
by | |66| , |67t and reviewed in [ |60| . But the key point can be understood simply: 
there is less structure in CHDM models on small scales because the growth 
rate of cold dark matter fluctuations is reduced on the scales where free 
streaming has wiped out neutrino fluctuations. Let us define the fluctuation 
growth rate / by 

fW . « . (9) 

a log a 

where S(k) is the amplitude of the fluctuations of wave number k — 2tt/ X in 
cold dark matter, and as usual a — 1/(1 + z) is the scale factor. For = 1 
CDM fluctuations, the growth rate f — 1- This is also true for fluctuations 
in CHDM, for k sufRciently small that free-streaming has not significantly 
decreased the amplitude of neutrino fluctuations. However, in the opposite 
limit k — >oo [||J60|, 



/oo = (v/l + 24a- l)/4 « 12,"-'^ , (10) 

assuming that f2c + Sli^ = 1. For example, for J?^ = 0.2, foo — 0.87. Even 
though the growth rate is only a little lower for these large-fc (i.e., short- 
wavelength) modes, the result is that their amplitude is decreased substan- 
tially compared to longer-wavelength modes. 

The next step in determining the implications for structure formation is 
to work out the effects on nonlinear scales using N-body simulations. This 



10 Primack and Gross 



is harder for Cold+Hot models than for CDM because the higher velocities 
of the neutrinos require more particles to adequately sample the neutrino 
phase space. The simulations must reflect the fact that the neutrinos initially 
have a redshifted Fermi-Dirac phase space distribution [Q. These CHDM 
simulations were compared with observational data using various statistics. 
CHDM with ill, — 0.3, the value indicated by approximate analyses | |6^ , |69{ , 
was shown to lead to groups of galaxies having substantially lower velocity 
dispersions than CDM, and in better agreement with observations But 
it also leads to a Void Probability Function (VPF) with more intermediate- 
sized voids than are observed |0 . This theory had so little small-scale power 
that a quasi-linear analysis using the Press-Schechter approximation showed 
that there would not be enough of the high-column-density hydrogen clouds 
at high redshift z ~ 3 known as damped Lyman-a systems |7^ , [73 74 . But 
CHDM with [2^, — 0.2 suppresses small-scale fluctuations less and therefore 
has a better chance of avoiding this problem j?^. Simulations |76) showed 
that this version of CHDM also has a VPF in good agreement with observa- 
tions ||77|| . The group velocity dispersions also remained sufficiently small to 
plausibly agree with observations, but it had become clear that the N-body 
simulations used lacked sufficient resolution to identify galaxies so that this 



statistic could be measured reliably |78 



A resolution problem also arose regarding the high-redshift damped Lyman- 
a systems. Earlier research had been based on the idea that these systems 
are rather large disk galaxies in massive halos |7^, but then high- resolution 
hydrodynamical simulations jSOt showed that relatively small gaseous pro- 
togalaxies moving in smaller halos provide a good match to the new, de- 
tailed kinematic data |Q. It thus appeared possible that CHDM models 
with Qi, ^ 0.2 might produce enough damped Lyman-a systems. With the 
low Hubble parameter h ~ 0.5 required for such S7„i — 1 models, the total 
neutrino mass would then be ^ 5 eV. 

While neutrino oscillation experiments can determine the differences of 
squared neutrino masses, as we will briefly review next, cosmology is sensitive 
to the actual values of the neutrino masses — for any that are larger than 
about 1 eV. In that case, cosmology can help to fill in the neutrino mass 
matrix. 

One example of this is the fact that if the hot DM mass is roughly evenly 
shared between two or three neutrino species, the neutrinos will be lighter 
than if the same mass were all in one species, so that the free streaming length 
will be longer. A consequence is that, for the same total neutrino mass and 
corresponding f2^, the power spectrum will be approximately 20% lower on 
the scale of galaxy clusters if the mass is shared between two neutrino species 
[|]. Since the amplitude and "tilt" n of the power spectrum in CDM-type 
models is usually fixed by comparison with COBE and cluster abundance, 
this has the further consequence that higher n (i.e., less tilt) is required when 
the neutrino mass is divided between comparable-mass neutrino species. Less 
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tilt means that there wiU be more power on small scales, which appeared to 
be favorable for the CHDM model, for example because it eased the problems 
with damped Lyman-a systems ■ 

5 Evidence for Neutrino Mass from Oscillations 

There is mounting astrophysical and laboratory data suggesting that neutri- 
nos oscillate from one species to another ||2^ , which can only happen if they 
have non-zero mass. Of these experimental results, the ones that are regarded 
as probably most secure are those concerning atmospheric neutrino oscilla- 
tions from Super-Kamiokande (see the chapter by John Learned) and solar 
neutrinos from several experiments (see the chapter by Wick Haxton). But 
the experimental results that are most relevant to neutrinos as hot dark mat- 
ter are from the Liquid Scintillator Neutrino Detector (LSND) experiment at 
Los Alamos (see the chapter by David Caldwell). 

Older Kamiokande data |Q showed that, for events attributable to at- 
mospheric neutrinos with visible energy £' > 1.3 GeV, the deficit of i'^ in- 
creases with zenith angle. The Super-Kamiokande detector has confirmed 
and extended the results of its smaller predecessor js^. These data imply 
that ^ Vt oscillations occur with a large mixing angle sin^ 2Q > 0.82 
and an oscillation length several times the height of the atmosphere, which 
implies that 5 x IQ-"* < Aml^ < 6 x lO^^ eV^ (90% CL). (Neutrino oscil- 
lation experiments measure not the masses, but rather the difference of the 
squared masses, of the oscillating species, here Arn^^ = \m{vr)'^ ~ "^(^m)^I-) 
This in turn implies that if other data requires either v^j_ or Vr to have large 
enough mass ( 1 eV) to be a hot dark matter particle, then they must be 
nearly equal in mass, i.e., the hot dark matter mass would be shared between 
these two neutrino species. Both the new Super-Kamiokande atmospheric z/g 
data and the lack of a deficit of 9^ in the CHOOZ reactor experiment 
make it quite unlikely that the atmospheric neutrino oscillation is —> Vg^. 
If the oscillation were instead to a sterile neutrino, the large mixing angle 
implies that this sterile species would become populated in the early universe 
and lead to too much ^He production during the Big Bang Nucleosynthesis 
epoch (Sterile neutrinos are discussed further below.) It may be pos- 
sible to verify that v^, —> Vr oscillations occur via a long-baseline neutrino 
oscillation experiment. The K2K experiment is looking for missing due to 
t'p ^ Vt oscillations with a beam of z^^ from the Japanese KEK accelerator 
directed at the Super-Kamiokande detector, with more powerful Fermilab- 
Soudan and CERN-Gran Sasso long-baseline experiments in preparation, the 
latter of which will look for r appearance. 

The observation by LSND of events that appear to represent v^J_ —>■ 
oscillations followed byPe+?'^"- + e~''i >I?-|-7, with coincident de- 

tection of e+ and the 2.2 MeV neutron-capture 7-ray, suggests that Avr?^^ > 
^7\ . The independent LSND data suggesting that z/^ — > Ve oscillations 
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are also occurring is consistent with, but has less statistical weight than, 
the LSND signal for i/^ — > oscillations. Comparison of the latter with 
exclusion plots from other experiments allows two discrete values of Z\m^g, 
around 10.5 and 5.5 eV^, or a range 2 eV^ ~ 0-2 eV^. The lower 

limit in turn implies a lower limit rrii, ^ 0.5 eV, or f2^ }t 0.01(0.65//i)^. This 
would imply that the contribution of hot dark matter to the cosmological 
density is at least as great as that of all the visible stars i7, ~ 0. 0045(0. 65//i) 
Such an important conclusion requires independent confirmation. The 
KArlsruhe Rutherford Medium Energy Neutrino (KARMEN) experiment has 
added shielding to decrease its background so that it can probe a similar re- 
gion of ^JTi^jg and neutrino mixing angle; the KARMEN results exclude a 
significant portion of the LSND parameter space, and the numbers quoted 
above take into account the current KARMEN limits. The Booster Neutrino 
Experiment (BOONE) at Fermilab should attain greater sensitivity. 

The observed deficit of solar electron neutrinos in three different types 
of experiments suggests that some of the i^e undergo Mikheyev-Smirnov- 
Wolfenstein matter-enhanced oscillations Ve — to another species of neu- 
trino Vx with Aml^ Ki 10~^ eV^ as they travel through the sun or pos- 
sibly "Just-So" vacuum oscillations with even smaller Am\^ [Q. The LSND 

—> Vf. signal with a much larger Aml^ is inconsistent with x = ^, and the 
Super-Kamiokande atmospheric neutrino oscillation data is inconsistent with 
X = T. Thus a fourth neutrino species Vs is required if all these neutrino oscil- 
lations are actually occurring. Since the neutral weak boson Z'^ decays only 
to three species of neutrinos, any additional neutrino species could not 
couple to the Z", and is called "sterile." This is perhaps distasteful, although 
many modern theories of particle physics beyond the standard model include 
the possibility of such sterile neutrinos. The resulting pattern of neutrino 
masses would have Ve and very light, and m{v^) « miyr) ~ (^™e/i)^^^i 
with the and Vt playing the role of the hot dark matter particles if their 
masses are high enough . This neutrino spectrum might also explain how 
heavy elements are synthesized in core-collapse supernova explosions Josf . 
Note that the required solar neutrino mixing angle is very small, unlike that 
required to explain the atmospheric z/^ deficit, so a sterile neutrino species 
would not be populated in the early universe and would not lead to too much 
^He production. 

Of course, if one or more of the indications of neutrino oscillations are 
wrong, then a sterile neutrino would not be needed and other patterns of 
neutrino masses are possible. But in any case the possibility remains of neu- 
trinos having large enough mass to be hot dark matter. Assuming that the 
Super-Kamiokande data on atmospheric neutrinos are really telling us that 
Ufj^ oscillates to i^r, the two simplest possibilities regarding neutrino masses 
are as follows: 
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A) Neutrino masses are hierarchical hke aU the other fermion masses, in- 
creasing with generation, as in see-saw models. Then the Super-Kamiokande 
Arri^ « 0.003 implies m(i^^) « 0.05 eV, corresponding to 

f]^ = 0.0013(TO^/0.05eV)(0.65//i)2 . (11) 

This is not big enough to affect galaxy formation significantly, but it is an- 
other puzzling cosmic coincidence that it is close to the contribution to the 
cosmic density from stars. 

B) The strong mixing between the mu and tau neutrinos implied by the 
Super-Kamiokande data suggests that these neutrinos are also nearly equal 
in mass, as in the Zee model and many modern models |^,^ (although 
such strong mixing can also be explained in the context of hierarchical models 
based on the SO(IO) Grand Unified Theory |9^). Then the above f}^, is just 
a lower limit. An upper limit is given by cosmological structure formation. 
In Cold -|- Hot Dark Matter (CHDM) models with Hm = 1, we saw in the 
previous section that if il^, is greater than about 0.2 the voids are too big 
and there is not enough early structure. In the next section we consider the 
upper limit on ili, if ilm ~ 0.4, which is favored by a great deal of data. 



6 Cold plus Hot Dark Matter and Structure 
Formation: firn ~ 0.4 

We have already mentioned that the ilm — 1 CHDM model with — 0.2 
was found to be the best fit to nearby galaxy data of all cosmological models 
But this didn't take into account the new high-z supernova data and 
analyses leading to the conclusion that I?a ~ ^matter ~ 0.2, nor the 
new high-redshift galaxy data. Concerning the latter, Somerville, Primack, 
and Faber found that none of the J7m — 1 models with a realistic power 
spectrum (e.g., CHDM, tilted CDM, or tCDM) makes anywhere near enough 
bright z ~ 3 galaxies. But we found that /ICDM with Sim ~ 0.4 makes about 
as many high-redshift galaxies as are observed |9^. This J7m value is also 
implied if clusters have the same baryon fraction as the universe as a whole: 
fim ~ ^b/ fb ~ 0-4, using for the cosmological density of ordinary matter 
f2b = 0.019/1^2 ra and foj. the cluster baryon fraction = 0.06ft.-3/2 pf 



from X-ray data or /f, = 0.077h ^ from Sunyaev-Zel'dovich data |100|. An 



analysis of the cluster abundance as a function of redshift based on X-ray 



temperature data also implies that fim « 0.44 ± 0.12 [101 102]. Thus most 
probably Sim is ^ 0.4 and there is a cosmological constant Ha ^ 0.6. In the 
1984 paper that helped launch CDM iQ, we actually considered two models 
in parallel, CDM with f2m = 1 and ylCDM with i7„ = 0.2 and f2A = 0.8, 
which we thought would bracket the possibilities. It looks like an /ICDM 
intermediate between these extremes may turn out to be the right mix. 

The success of — 1 CHDM in fitting the CMB and galaxy distribu- 
tion data suggests that fiat low-J?^ cosmologies with a little hot dark matter 
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Power Spectrum for n=l ACDM and ACHDM 




10-3 iQ-s 10-1 IQO IQl 

k (h Mpc-i) 



Fig. 1. Nonlinear dark matter power spectrum vs. wavenumber for ylCDM and 
ylCHDM models with O^/n^ = 0.05,0.1,0.2,0.3. Here = 0.4, the Hubble 
parameter h — 0.65, there is no tilt (i.e., n = 1), and the bias b — 0.85. Note tha t 
in this and the next Figure we "nonlinearized" all the model power spectra |l03t , 
to allow them all to be compared to the APM data (the small "wiggles" in the 
high-/2j, power spectra are an artifact of the nonlinearization procedure). 
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Fig. 2. Nonlinear dark matter power spectrum vs. wavenumber for 12 ylCHDM 
models with Ni, = 2 massive neutrino species and Hubble parameter h = 0.65, with 
tilt and erg determined by COBE + ENACS cluster abundance. The bias chosen 
for these models is that which minimizes over the entire range of available APM 
data. 
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Fig. 3. CMB anisotropy power spectrum vs. angular wave number for the same 
models as in Figur e 2. Tlie data pl otted are from COBE and three recent small- 
angle experiments [ p^4|Jl05| , p^ , p^ |. 
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be investigated in more detail. We have used CMBFAST IQ to examine 
ylCHDM models with various h, and flj,, assuming i7f, = 0.019/i^^. Fig- 
ure 1 shows the power spectrum P{k) for ylCDM and a sequence of ylCHDM 
models with increasing amounts of hot dark matter, compared to the power 



spectrum from APM |108 . Here we have fixed f?™ = 0.4 and Hubble param- 
eter h = 0.65. All of these models have no tilt and the same bias parameter, 
to make it easier to compare them with each other. As expected, the large- 
scale power spectrum is the same for all these models, but the amount of 
small-scale power decreases as the amount of hot dark matter increases. 

In Figures 2 and 3 we consider a sequence of twelve ylCDM and ylCHDM 
models with h = 0.65, fim = 0.3,0.4,0.5,0.6, and f?^//2,„ = 0,0.1, and 0.2. 
We have adjusted the amplitude and tilt n of the primordial power spectrum 
for each model in order to match the 4 -yea r CO BE a mplitude and the ENACS 
differential mass function of clusters |10£l (cf. |11C|). (We checked the CMB- 
FAST calculation of ylCHDM models against Holtzman's code used in our 
earlier investigation of ylCH DM mode ls [0 . Our results are also compatible 
with those of recent studies | 111, 11 2 1 in which n = 1 models were consid- 
ered. But we find that some ylCDM and ylCHDM models require n > 1, 
called "anti-t ilt" , and it is easy to create cosmic inflation models that give 
n > 1 — cf. |113|.) In all the ylCHDM models the neutrino mass is shared 
between N^, = 2 equal-mass species — as explained above, this is required by 
the atmospheric neutrino oscillation data if neutrinos are massive enough to 
be cosmologically significant hot dark matter. (This results in slightly more 
small-scale power compared to N^, = 1 massive species, as explained above, 
but the Ni, = 1 curves are very similar to those shown.) In Ref. |114| we have 
shown similar results for Hubble parameter h = 0.6, and also plotted the best 
CHDM and ylCDM models from |^ . Note that all these Figures are easier to 
read in color; see the version of this paper on the Los Alamos archive [ [115| . 

Of the ylCHDM models shown, for Sim = 0.4 — 0.6 the best simultaneous 
fits to the small-angle CMB and the APM galaxy power spectrum data [108| 
are obtained for the model with Qu/ Sim = 0.1, and correspondingly ■m[Vfj) k, 
m{i'r) ~ 0.8 — 1.2 eV for h = 0.65. For Sim < 0.4, smaller or vanishing 
neutrino mass appears to be favored. Note that the anti-tilt permits some of 
the ylCHDM models to give a reasonably good fit to the COBE plus small- 
angle CMB data. Thus, adding a little hot dark matter to the moderate- 
ylCDM models may perhaps improve somewhat their simultaneous fit to the 
CMB and galaxy data, but the improvement is not nearly as dramatic as was 
the case for Sim = 1- 

It is apparent that the ylCDM models with Sim = 0.4, 0.5 have too much 
power at small scales (fc ^ lh~^ Mpc), as is well known |117,118| — although 
recent work ]11S[ | suggests that the distribution of dark matter halos in the 
Sim = 0.3, h = 0.7 ylCDM model may agree well with the APM data. On the 
other hand, the ylCHDM models may have too little power on small scales 
— high-resolution ylCHDM simulations and semi-analytic models of early 
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— I 1 , , 1 r- 

(a) N=l (species) 




Fig. 4. Constraints on tlie neutrino mass assuming (a) N^, = 1 massive neutrino 
species and (b) Ni, = 2 equal-mass neutrino species. Tlie lieavier wei ght curves 
show the effect of including the Lyman-alpha forest constraint. (From [116|; used 
by permission.) 



galaxy formation may be able to clarify this. Such simulations should also 
be compared to data from the massive new galaxy redshift surveys 2dF and 
SDSS using shape statistics, which have been shown to be able to discriminate 
between CDM and CHDM models |l20 |. 

Note that all the ylCDM and ylCHDM models that are normalized to 
COBE and have tilt compatible with the cluster abundance are a poor fit 
to the APM power spectrum near the peak. The ylCHDM models all have 
the peak in their linear power spectrum P{k) higher and at lower k than the 
currently available data (e.g., from APM). Thus the viability of ylCDM or 
ylCHDM models with a power-law primordial fluctuation spectrum (i.e., just 
tilt n) depends on this data/analysis being wrong. In fact, it has recently 
been argued [121| that because of correlations, the errorbars underestimate 
the true errors in P{k) for small k by at least a factor of 2. The new large- 
scale surveys 2dF and SDSS will be crucial in giving the first really reliable 
data on this, perhaps as early as next year. 
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The best published constraint on in ylCHDM models is ]116[ | . Figure 4 
shows the result of their analysis, which uses the COBE and cluster data much 
as we did above, the P{k) data only for 0.025(/i/Mpc) < fc < 0.25(/i/Mpc), 
the constraint that the age of the universe is at least 13.2±2.9 Gyr (95% C.L.) 
from globular clusters [122|, and also the power spectrum at high redshift 
z ^ 2.5 determined from Lyman-a forest data. The conclusion is that the 
total neutrino mass rrii, is less than about 5.5 eV for all values of f2m, and 
rrii, ^ 2.4(i7,„/0.17 — 1) eV for the observationally favored range 0.2 < < 
0.5 (both at 95% C.L.). Analysis of additional Lyman-a forest data can allow 
detection of the signature of massive neutrinos even if rrii, is only a fraction 
of an eV. Useful constraints on f2^ will also come from large-scale weak 
gravitational lensing data |122] combined with cosmic microwave background 
anisotropy data. 

JRP acknowledges support from NASA and NSF grants at UCSC. 
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